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' Abstract 

CN 

We study a relation recently derived by K. Kondo at zero momentum between the 
D ' Zwanziger's horizon function, the ghost dressing function and Kugo's functions u and w. 

, We agree with this result as far as bare quantities are considered. However, assuming the 

validity of the horizon gap equation, we argue that the solution w{0) = is not acceptable 
since it would lead to a vanishing renormalised ghost dressing function. On the contrary, 
when the cut-ofF goes to infinity, u{0) oo, w{0) — oo such that u(0) + w(0) — 1. 
' (— I 'I Furthermore w and u are not multiplicatively renormalisable. Relaxing the gap equation 

pi i' allows w{0) = with u{0) ^ — 1. In both cases the bare ghost dressing function, F{0, A), 

JL I goes logarithmically to infinity at infinite cut-ofF. We show that, although the lattice 

results provide bare results not so different from the F{0,A) = 3 solution, this is an 
(~| I accident due to the fact that the lattice cut-offs lie in the range 1-3 GeV~^. We show that 

the renormalised ghost dressing function should be finite and non-zero at zero momentum 
and can be reliably estimated on the lattice up to powers of the lattice spacing ; from 
^ ' published data on a 80^ lattice at /3 = 5.7 we obtain Fj^{0, i^l = 1.5 GeV)~ 2.2. 

in 
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^ ■ 1 Introduction 



X 



Kondo has derived in recent papers [1,2] a relation between the fc = values of the 
ghost dressing function F{k), Zwanziger's horizon function h{k), Kugo's function u{k), and 
' an additional function w{k). Applying to this relation Zwanziger's horizon gap equation 

and assuming that w{0) = he derives the surprising result that n(0) = —2/3. This is 
surprising as so simple constraints on bare quantities are rare. We know only the case of the 
electric charge which benefits of the Ward identity. This surprising result deserves some closer 
investigation, even more so as lattice results are not so far from it as we shall see. Indeed it 
has given rise to several pubhcations and there is far from a consensus on this matter [3-6]. 

To understand better the issue we try in this note to reconsider every point of the dis- 
cussion from first principles. Our starting point is a set of relations between the functions 
we have just mentionned. They concern bare quantities, which imposes to use a finite ultra- 
violet cut-off, else we would have to deal with divergent quantities. In section [5] we propose 
a faster derivation of these relations. If one assumes the validity of the Zwanziger horizon 
gap equation this boils down to very simple relations giving u{Q) and w{0) as functions of 
F(0). We then discuss whether u(0) can be -2/3 or any finite quantity. We argue that it is 
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not possible if we assume F to be multiplicatively renormalisable, which nobody would deny. 
In section [3] we use lattice QCD and ghost propagator Dyson-Swinger equation (GPDSeq) to 
get numbers. From the GPDSeq at small momentum we find that ratio of the the bare (resp. 
renormalised) ghost dressing funtions at small and zero momentum, assuming the latter to be 
finite, is essentially cut-off and renormalisation point independent. We extract an estimate of 
the renormalised -F/j(0). Finally we discuss the status of the Zwanziger horizon gap equation 
on the lattice. Convinced that it has no reason to be valid, we generalise the result of section 
[2] for a more general case. 

2 Ghost dressing function, horizon function, u and w 

The discussion which follows deals with bare quantities. These are singular and need a 
regulator, or cut-off, which we will call A (in the lattice case, this regulator is a~^, a being 
the lattice spacing). The dependence in A will often be kept implicit, to avoid heavy notations, 
but is always understood speaking of bare quantities. Renormalised quantities will be marked 
by the index R. There is no need to specify the renormalisation scheme being used, since 
our results do not depend on a particular choice; however, regarding lattice results, we shall 
refer as usual to the MOM scheme. 

2.1 Gribov-Zwanziger action 

In [1,2], it has been claimed that three-point and four-point functions for gluon and ghost 
fields can be related in such a manner that the Zwanziger horizon condition strongly constrains 
the ghost propagator and the ghost-gluon vertex. 

It is well known, since Gribov's famous paper [8], that the gauge fixing procedure in QCD 
using the standard Faddeev-Popov procedure is not unambiguous. It leads to a discrete 
set of solutions, named "Gribov copies". One solution, proposed by Zwanziger [9] , which 
aims at restricting the Gribov copies within the Gribov Horizon, consists in using the Gribov- 
Zwanziger partition function in Landau gauge, 

Z-y = J [DA]6{dA) det(M) e-SYM+ 7/^°^^^) ^ (1) 

for the D-dimensional Euclidean Yang-Mills theory, where 5ym stands for the Yang-Mills 
action, M is the Faddeev-Popov operator, 

= -d^Df = -d, {d,5-'> + gr^'^A'^^ (2) 

and h{x) is the Zwanziger horizon function, 

h{x) = I d^y gr^'^Al{x){M-^r{x,y)grf-Al{y) ; (3) 

that restricts the integration over the gauge group to the first Gribov region, provided that 
the Gribov parameter, 7, is a positive number that is to be determined by solving the so-called 
gap equation: 

{h{x)), = {N^-l) D. (4) 

The horizon function is a bare quantity depending on the cut-off parameter, as 7 also does 
through the implementation of the horizon condition that requires that the gap equation be 
solved for every cut-off value. We will postpone the discussion of this gap equation to a later 
section. In any case, the horizon function is a well defined bare quantity and it is relevant to 
first derive, independently of eq. ^ the relation its v.e.v. has with other quantities. 
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2.2 Relating h{0),u{0),w{0) and F(0) 

In this subsection we propose a simplified derivation of Kondo's relation (cf [1, 2]) which 
relates the v.e.v. of the horizon function h{0) to the ghost dressing function at vanishing 
momentum and the Kugo-Ojima parameters. Then, contrarily to Kondo, we will add no 
assumption about the Kugo-Ojima parameters but simply combine Kondo's relation with 
the one discovered by Kugo between the ghost dressing in Landau gauge with these Kugo- 
Ojima parameters and discuss about their general implications. 

{hmk=o = Hm I d^x I d^'y {gr''Al{x){M~^r{x,y)grf'Al{y)) e'^^^-y^ 
= lim:^ j d^x j d^y { (^gr^^A^ c^) (^gr^^Af^^^ ^eM^-y) 
= lim f d''{x-y){(gr^'^AlA (gr^'^AlA )e^''<^-y^ 
= { [gr^'^A';, c j [gr^'Al 7f) ),.^, (5) 
where we use the simplified notation: 

{{...){...)), ^ j d^{x - y){{. . . U. . . )y) e''<^-y^ (6) 

that was introduced in ref. [1,2] and that will be followed from now on. To establish eq. ([5]), 
nothing is needed but the relation between the inverse Faddeev-Popov operator and the ghost 
and anti-ghost fields and the translational invariance. Define then the function ti(fc^), the 
value of which at vanishing momentum gives the Kugo-Ojima parameter, as 

( {Dfc') {gr'^Atlf) )k =- ^IJ""' u(k^) ; (7) 

where k'^Sj^^{k) = k'^S^y — k^ku and the transversality is guaranteed by the well-known iden- 
tity: 

( (S^^c") [gf^'^At-/) )k = -ik,{ [Dfc^) [gr^'Atc') )fc = . (8) 

Now, by merely invoking the definitions of u (eq. ([7|)) and of the covariant derivative, 

Df^5-'d, + gr'A'^^, (9) 
acting on the ghost and anti-ghost fields, one obtains 

{(gr''Alc')(gf'fAl-cf))u = {{Dl'c'){gf'fAl-cf))k- { d^d^ [gf'f Alcf)),, m 

-SlAkW'uik^) -^k,{c'^{gf^fAlcf)), 
which proves that the transversal part of the l.h.s. of eq. (jlU|) is given by eq. ([7]), 

C' {{aP'^A'^.c^) {gf'fAlcf)), = -6lAk)5'^Mk^) . (11) 
while the longitudinal part can be written as follows: 

k,{(gr''Alc')(gf'fAl-cf))k = ik\c'^{gf'fAl-cf))u 

= tk^{^^{e' {gf'fAl-cf))^; (12) 
5'^'^'F(A;2) 
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where F{k'^) is the bare ghost propagator dressing function and IPI notes the one-particle 
irreducible contribution to the v.e.v. obtained through the amputation of the external ghost 
leg. Let us then define the function w{k'^) in order to parametrize this longitudinal contribu- 
tion to eq. through 

(c« [gf^'^lcf) )r = Hk") + wi^)) . (13) 

This definition is equivalent to the one given in terms of diagrams in ref. [1,2]. Note also that 
w{0) was taken to be in the seminal work by Kugo and Ojima. 
Taking together eqs. ([TO]) and (fT3]) one gets: 

{{gr'''Kc^){9f''^Al-cf)), =-5-'^[5lMk') + F{k^)^-^ {u{k^)+w{k^))^ , (14) 
and for the v.e.v of the horizon function [1,2]: 



{hmk=o 1 



= [{D - l)u(O) + F(0) {u{Q) + wm ] . (15) 

From now on we make explicit the dependence on the cut-off, A, of all the bare quantities, 
which generally will diverge in the infinite cut-off limilH. 

Kugo has shown in refs. [7, 10] that the Landau gauge condition, d^A^^^ = 0, can be 
exploited to give: 

(l + w(0,A) +u;(0,A))F(0,A) = 1 . (16) 

This result can also be easily derived from the ghost-propagator Dyson-Schwinger equation 
which, in Landau gauge, can be written as: 



safe 



= l + u{k\K)+w{k\K) (17) 

Then, the two equations ([15]) and p6]) can be combined to obtain, without any hypothesis 
about u and w. 



u{0,A) 
u;(0,A) 



F(0,A)-1 D 
D-1 D - 

-l-n(0,A) + 



{hmk=o 
D{N^ - 1)_ 



F(0,A) 
F(0,A) + (I)-2) 



(18) 



D-1 



+ 



+ 



D 



F(0,A) D-l 



k=0 



D{N^ - 1) 



as general solutions for the Kugo-Ojima parameters, n(0, A), and w{0,A) in terms of the bare 
ghost dressing function at vanishing momentum. If, in addition, we assume that the gap 
equation eq. ^ is satisfied the square bracket in eq. (fT8|) is equal to 1, independently of the 
cut-off and we get 



^It is well-known that the ghost dressing function diverges logarithmically at infinite cut-off in the UV 
momentum domain. 
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u(0,A) 



F(0,A) -1-D 
D-1 



(19) 



u;(0,A) 



l-u(0,A) + 



1 



F(0,A) -2 
D-1 



+ 



1 



In fig. [H the solutions of eq. (fT9]l are plotted as functions of F(0, A). It is obvious from 
eq. (fT^ (see also fig. [T]) that, had we required w{0,A) = 0, the solution proposed in ref [1,2] 
would emerge: u(0. A) = —2/3 and -F(0, A) = 3, for = 4. However we shall present in the 
next subsection the arguments which lead us to believe that no solution implying a cut-off 
independent and finite valu^ for F{0) can be accepted. 

2.3 Constraints from renormalisability 

In this section we assume the validity of the relation ^ Let us start from the basic equation 



where Fji is the renormalised dressing function in the infinite cut-off limit and n some 
positive number. For any quantity which is multiplicatively renormalisable, a field-theory 
non-perturbative renormalisation scheme (in particular, those applied in lattice field theory) 
implies a relation of this kind where the crucial point is that the cut-off dependence is an in- 
verse power of the cut-off and cannot behave like some power of the cut-off's logarithm [13,14]. 
Now, we know from perturbation theory how Z^{fi'^,A) depends on the cut-off; choosing a 
fixed value Aq and sending A ^ oo, diverges logarithmically in the infinite cut-off limit, 
regardless of the renormalisation procedure : 



Although this behaviour is quite general, the specific value ^ of the exponent is valid only 
in the case A'^ = 3, Nf = 0. This provides us with two main objections for a finite bare ghost 
dressing function: 

• A finite value of F(0, A) requires, through eq. (f20]) . that F/j(0,/x^) be zero; taking into 
account the fact that the subdominant terms, which vanish as A — > oo, are supposed 
not to be logarithmic contributions but, at least, of the order of 1/A we are forced to 
conclude that zero is the only allowed finite value that the bare ghost dressing function 
can hit for eq. ()20p to be consistently satisfied. 

• We can apply eq. (|29ap which will be discussed later and implies at small a cut-off 
independent factor, decreasing with q^, which multiplies -F(0, A) and as well Fr(0,//^). 
Recalling that the path integration has been limited by the Zwanziger procedure to a a 
region in which the Faddeev-Popov operator is positive , we see that if Fh(0,/x2) = 0, 
Fr can only assume the value throughout some range of g^, which sounds weird. Even 
more, there are numerical evidences that the ghost dressing function F{q'^, A) decreases 
for all q^. Then Fji{q'^,^'^) = should hold for any q^. 

Therefore, the only way out we see is that the bare ghost dressing function diverges 
logarithmically in the infinite cut-off limit and that multiplication by Zg provides a strictly 
positive renormalised value. Then, eq. ()19p can be rewritten as: 

^This finite value is independent of the number of colors and, provided that w{0, A) = 0, should be the 
same whichever regularized Yang-Mills action we use (any lattice action, for instance) or even including any 
non-zero number of quark fiavours for the action. 




(20) 




9/44 




(21) 
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Figure 1: The solutions for m(0,A) and w(0, A) given by eq. plotted as a function of F(0,A). 
This plot can be understood as a function of Z3 for a given non zero value of Fr(0, /i^). Then the 
infinite cut-off limit is the limit at infinity of the horizontal axis. The particular solution proposed in 
ref. [1,2] (black circles), obtained by imposing w{0, A) = 0, corresponds to the intersection of u + w and 
u. It cannot hold when 00. The current lattice solutions for the bare ghost dressing functions 

at vanishing momentum lie inside the green dotted square (see fig. [3]). This apparent approximate 
agreement is misleading and due to the moderate cut-off value on the lattices. 



Fn{0,^l') = Z^'{f,',A){{D-l)u{0,A)+D + l 



u{0,A) +w{0,A) 



(D-l) 



-1 + 



Z^\fi\A)u{0,A) + O 



log A 



(22) 



log A 



An important consequence of the first of equations ()22|) is that the function u{0, A) cannot be 
multiplicatively renormalised. The fact that, when multiplied by Z^^ , it gives a finite result 
in the infinite cut-off limit does not suffice. As we have already recalled, it has to differ from 
its asymptotic value by inverse powers of the cut-off, which is obviously wrong in eqs. ()22p 
where Z^ n(0. A) only converges up to inverse powers of the logarithm of the cut-off. The 
same is true for w^O). This is not surprising since they are defined by the insertion of the 
composite operators shown in eqs. (|13|7|) . 

Only the combination 1 -|- u(0, A) -|- w(0, A) vanishes logarithmically as A —> 00 so that. 



Z3ifi^A)il + uiO,A) + wiO,A)) 



1 



i^/?(0,/.2; 



+ 



(23) 



while both u{0,A) and w{0,A) diverge but their divergences cancel in eq. (f23l) . Thus, as done 
in [2, 11], one can consider 1 + u + w to he renormaliseci^ by Z3. However, let us repeat, 
1 + u and w cannot be separately renormalised and 7^(0, A) = cannot be accepted since 
in conjunction with equations ([!]) and (jlSp it provides a finite u{0,A) = —2/3 and a finite 
-F(0, A) = 3 which we have shown above to be forbidden. 



* Renormalised in the sense that the (logarithmically vanishing) cut-off dependence can be killed at a given 
renormalisation point up to vanishing powers of the cut-off 



6 



3 Collecting and extrapolating bare ghost lattice data 



Lattice simulations first provide us with estimates for bare quantities (correlation func- 
tions) in the lattice regularization scheme, where the role of the regularization cut-off is played 
by the inverse of the lattice spacing, . In present simulations is moderate, ranging 
from ~ 1 GeV~^ for /3 = 5.7 up to ~ 3.5 GeV~^ for /3 = 6.4. Those bare quantities should be 
further renormalised by applying MOM-like schemes. That this multiplicative renormalisa- 
tion procedure works has been proven by Reisz in ref. [13]; the remaining corrections due to 
finite spacing (vanishing in the continuum limit) are considered to behave as powers of the 
lattice spacing. Those renormalised quantities are usually the reliable result of the simula- 
tions and the ones directly connected with physical predictions. On the contrary, the recent 
work [1,2] we discussed above supplies a prediction for a bare quantity: the bare ghost prop- 
agator dressing function. Therefore, the non-renormalised lattice estimates for this dressing 
function deserve by themselves a particular interest, as far as they could allow us to test that 
prediction. 

In the last few years, many works have been devoted, at least partially, to the lattice 
computation of the ghost propagator. They mainly follow ref. [15] in writing the Faddev- 
Popov operator lattice divergence: 

MiU) = -^V-DiU) (24) 

where the operator D reads 

Df,{U)ri{x) = ^ (U^{x)r,{x + /i) - 7]{x)U^{x) + r,{x + ft)Ul{x) - C/^(x)r/(x)) (25) 

Those definitions, complemented with conversion routines between the Lie algebra and the 
Lie group, allow for a very efficient lattice implementation. Some details about the procedure 
for the inversion of the Faddeev- Popov operator and some results will be found in [16]. 

The gauge fixing, in particular for Landau gauge, is a more delicate issue. A minimization 
of the functional 

Fu[g\='^eY,Y. (l- ^5(a;)t/^(x)5t(x + A)) (26) 

can be achieved by the use of some algorithm driving the gauge configuration to a local 
minimum of Fu\g\. The gauge configurations obtained in this way will lie in the first Gribov 
region but, in general, they do not reach the fundamental modular region defined as the set of 
absolute minima of -Ft/ [5] on all gauge orbits. A "best-copy" algorithm (basically consisting in 
choosing the gauge configuration providing the lowest minimum after several minimizations) 
has also been used as well as a procedure that essentially consists in a simulated annealing 
technique and is claimed to reach gauge-functional values closer to the global minima than 
the standard approach (see for instance [17,18] and references therein). Figure [3] presents 
together results collected from ref. [17] (for very big lattice-volume simulations with the 
simulated annealing gauge-fixing) and data from ref. [16,19,20] (obtained using the standard 
gauge-fixing); it shows only a weak dependence in the cut-off a~^. This is not surprising since 
one knows from eq. ()2ip that it should behave at leading log as which gives no larger 

effect than 2.5 % on the whole range of /?'s. 

In order to compare the lattice data with the previous results we need to extrapolate them 
to zero momentum. To carry this task out we derive now from the bare ghost propagator 
Dyson-Schwinger equation (GPDSeq) a small-momentum formula the coefficients of which 
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Figure 2: Bare zero- momentum gluon propagator estimated from different lattice data sets plotted 
in terms of the inverse of the lattice size in physical units. The data for the two bigger lattice volumes 
are taken from ref. [17], the smaller volume at /3 = 5.7 corresponds to ref. [26] and the others to 
refs. [16,20]. A hnear fit for the three data at /3 = 5.7 to extrapolate at infinite volume is presented 
as a solid red line. 

are fitted against the lattice data. The bare GPDSeq can be regularised and evaluated with 
the help of a subtraction procedure at two different momenta p and k, 

as explained in ref. [20,210. In this equation N is the number of colours, g{A) the bare, cut- 
off dependent, coupling constant and G stands for the gluon propagator dressing function 
and Hi is one of the form factors for the bare ghost-gluon vertex, 

ff'^i-q, k-q-k)= -igor"' ( qMq, k) + (q - k),H2{q, k) ) , (28) 

that should be finite and only weakly dependent on the momenta by virtue of Taylor's non- 
renormalisation theorem [22] and that, consequently, is usually assumed to be constant with 
respect to the momenta. Such a bare (and cut-off dependent) GPDSeq can be numerically 
solved, as was done in ref. [20], with the help of the lattice gluon propagator estimate to 
be inserted in the integral in eq. p7p . It is known (cf. [20]) that the solutions can belong 
to 2 different types : while the generic solution goes to a finite non-zero limit in the infra- 
red there exists also an exceptional one (for a given value of the coupling) which diverges 
as XjyfW' in this same limit. The lattice simulations clearly favour the first type and it is 
implied in the coming discussion that we are in this situation. The only unknown ingredient 
is the (constant) value for the ghost-gluon vertex form factor, ffi(A). In addition, following 
ref. [21], one can derive a small-momentum expansion for its solution 



In these references we dealt with renormalised quantities but everything applies straightforwardly for bare 
ones. 
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k (GeV) 



Figure 3: Bare ghost dressing function estimated from different lattice data sets. The data for the 
two larger lattice volumes are taken from ref. [17] and the others from refs. [16,20]. the solid line is 
for the best fit with the small-momentum expansion in eq. ([29a|) with R{(3 = 5.7(80^)) and the dashed 
one stands for the best fit with i?(/3 = 5.7, oo), both computed as explained in the text. 



with: R{A) 



F(0,A) 



1 + 



NHi{A)R{A) 2 



, F(0,Af ( lim 
47r \k^o 



G{k'^,A) 



fc2 



lim 



A;2 



(29a) 
(29b) 



where ot is the coupling constant defined in the Taylor scheme [23] . Before turning to exploit 
this expansion we shall comment briefly on the various parameters it involves. 

Insofar as the gluon propagator reaches a finite non-zero value at vanishing momentum 
(as it appears to be true on the lattice), R{A) takes a finite value which depends on the 
cut-off by inverse powers. Thus, the dominant (in relative terms) q^.^igpgndent part of F 
in the vicinity of (i.e. the second term in the bracket of eq. (I29a|) ) does not require any 
renormalisation, 



A^oo F(0,A) 



l + NHi 



levr 



■\og{q^)+0{q^ 



(30) 



This (quasi-)independence of the slope with respect to the cut-off is of course important to 
ensure multiplicative renormalisability, since the latter demands that the g'^-dependence of the 
bare and renormalised Green's functions be the same up to negative powers of the cutoff. On 
the contrary the global multiplicative factor F(0, A) is known to be logarithmically divergent 
with A, which implies that its variation could be appreciable. 

Since R{K) is to be evaluated from lattice estimates of zero-momentum gluon and ghost 
propagators, one expects it to be sensitive to finite-volume artefacts, to which much attention 
should therefore be paid. In fig. [2] one notices the strong dependence of the zero-momentum 
gluon propagator on the lattice size, which implies an equally strong dependence for i?(A). 
When evaluated from the zero-momentum estimates for a 80^-lattice at /5 = 5.7 in ref. [17], 
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R{(3 = 5.7) takes on the value of 10.3 (while, for instance, from the data for a 32^-lattice 
at /3 = 5.8 in ref. [20], one would obtain R{(3 = 5.8) ~ 19). The bare vertex form factor 
(supposed to be constant) was indirectly estimated in ref. [20] for a 32'^-lattice at /? = 5.8 
and appeared to be Hi((3 = 5.8) ~ 1.2. 

Provided that finite-size and lattice spacing artefacts can be neglected for the bare ghost- 
gluon vertex, the values of R and Hi we have just determined can be used to attempt to 
describe the ghost dressing function at small momenta estimated from the 80^ lattice at 
f3 = 5.7 in ref. [17]. Then the only parameter in eq. (I29ap which remains to be determined 
by the best fit is the zero-momentum ghost dressing function. Actually, since a value for 
-F(0, A) is required in computing R{A) one has to proceed by iterations: the known value of 
G(^-'(0, A) and an initial guess of F{0,A) are inserted in eq. ()29bp to produce a first estimate 
of R{A). The latter is then used in eq. ()29ap to perform a new fit of the 80^-lattice data 
deprived of the few momenta with p < Att/L. It appears actually that, for all the lattice 
data sets plotted in fig. [3l only the momenta satisfying his condition are affected by sizeable 
lattice volume artefacts. This produces a new estimate of -F(O) and the process is iterated 
until it eventually converges to 

i2(/3 = 5.7, 80) = 10.3, F(/3 = 5.7, 80) = 3.50. (31) 

The fit is presented as a solid line in fig. El 

The impact of the finite-volume effect due to the lattice determination of R{A) can be 
approximatively estimated in the following way. First, the zero- momentum gluon propagator 
data for the three different lattice volumes at (3 = 5.7 in fig. [2] can be extrapolated up to 
infinite volume (we work only with data for the same (3, in order avoid any mixing between 
lattice-spacing and volume effects). This is the starting point to repeat the iterative proce- 
dure explained above and one gets R{(3 = 5.7, cxd) = 8.2 and F((3 = 5.7, oo) = 3.40. The 
extrapolation is also shown (with a dashed line) in fig. [3l 

As for the bare Kugo-Ojima parameter u(0. A), it is estimated to be of the order of 
— (0.6 — 0.8) in ref. [24,25] and very recently in ref. [11] by using a mixed approach, analogous 
to the one previously applied to solve eq. ((271) . in which DS equations are solved with the 
input of a lattice estimate of the gluon propagator. 

3.1 Horizon gap equation and lattice QCD 

A few words are in order to compare the two approaches we have considered in this note. 

The Gribov-Zwanziger (GZ) approach eq. ([T]) proposes a modification of the standard 
QCD action in order to limit the domain of the path integration to the domain within the 
Gribov horizon, in which the Faddeev-Popov (FP) operator is positive, i.e. all its eigenvalues 
are positive. 

Lattice QCD uses the genuine QCD action. Some algorithm minimises the functional 
which discretises the functional / d'^xAaAa- These algorithms differ, but they all stop at a 
local minimum. Local minimum means that all the second order derivatives are positive, i.e. 
that the FP operator is positive. 

Therefore the two approaches share the property that they limit themselves to the in- 
terior of the Gribov horizon. None of them manages to find the absolute minimum of that 
functional, i.e. to stay within the fondamental domain. There is also a "thermodynamic" 
argument claiming that one stays close to Gribov's horizon which means in a domain where 
the eigenvalues of the FP operator should be small. This also seems to be valid for both 
approaches. 

Now come the differences. The GZ approach gives some weight to the different Gribov 
copies within the Gribov horizon. The lattice algorithms give another one, which, moreover, 
presumably depends on the specific features of each algorithm: it has been argued that 
stimulated annealing [17, 18] leads to smaller values of the minimised functional. 
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— K(A) 

— u(0,A) 
• Kondo's solution 

(n(0,A) 




Figure 4: The same plot shown in fig. [T] but the gap equation is given here by eq. ([5^ . the sohd blue 
line being for the new factor k(A) in that equation, and w(0, A) is required to be zero, as explained 
in the text. Again, current lattice estimates lie inside the green dotted square. 



On the lattice it is possible in principle to compute the v.e.v. of the horizon function 
(/i(0)). Nothing imposes that it should be independent on the cut-off and we do not see any 
reason why it should verify the horizon gap equation eq. This is precisely a place where 
the differences we just mentionned could be visible. 

In order to discuss this situation let us define a factor k(A) such that 



(/i(0))fc=o 
The value k(A) = 

n(0,A) 
w{0,A) 



lim 



d^x {h{x))e 



Ik-x 



k(A) (iV^ -1) D 



(32) 



1 corresponds to eq. Eq. (fTSl) now reads 
F(0,A)-1 Dk{A) 



D-1 
-1 - n(0,A) 

F{0,A) + D 



+ 



D-l 
1 



F(0,A) 
2 



D - 1 



+ 



1 



+ 



Dk{A) 
D-l 



(33) 



F(0,A) 

If the gap equation eq. (H]) is relaxed in this way it becomes possible to have w{0) = 
as shown in [11, 12] in the Landau background gauge and assumed in [1,2]. The solution 
becomes 

D-2 D-l 



n(0,A) 



i^(0,A) 



D 
1 



+ 



D DF{0,A) 



(34) 



i^(0,A) ■ 

implying that k{A) oo and n(0, A) — > —1 when A ^ oo (see fig. H]). This does not change 
our major conclusion that no finite value of F{0, A) independent of A is acceptable. In 
particular it happens that u(0, A — > oo) converges to —1, nevertheless Kugo-Ojima's 
condition only emerges at the infinite cut-off limit and thus < Fr(0, /z^) < oo. 



Lattice measurements, which correspond to A ~ 1 
which, with u ~ (.6 — .8) and F(0)~ 3.5, leads to k(A) 
tell how k(A) depends on A. 



4 GeV, give an estimate of k(A) 
: 1.2 — 1.3. However this does not 
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4 Conclusions 



We have generalized the solution recently proposed by Kondo for the zero-momentum ghost 
dressing function and the Kugo-Ojima parameter, m(0, A), by deriving eqs. (jlSp where both 
Kugo-Ojima parameters, u{0,A) and w{0,A) appear written in terms of F{0,A) and the 
horizon function, {h{0,A)), at any finite cut-off. In particular, we have shown that the one 
relating u{0,A) and F(0,A), after applying the gap equation, eq. (jl]), is close to be verified 
by lattice estimates for A = ~ IGeV, but that this is a pure coincidence due to the small 
cut-off in lattice calculations. We have argued that neither u nor w can be multiplicatively 
renormalised. Indeed, from the anomalous dimension of the ghost propagator renormalisation 
constant we conclude that no solution with a cut-off independent bare F{0, A) is possible. 
If the gap equation is valid, for A — > oo, then n(0, A) — > oo and w{0,A) — oo such that 
n(0. A) -|- w{0, A) — > —1. If one relaxes the gap equation eq. (f32]l . one can satisfy wi^O, A) = 
with fi:(A) — > oo and n(0, A) — > —1. The Kugo-Ojima condition n(0) = —1 is asymptotically 
fulfilled for the bare u while the renormalised ghost dressing function is finite and non zero. We 
have argued that lattice QCD, notwithstanding some similarity with the Gribov-Zwanziger 
approach, has no reason to fulfill Zwanziger's horizon gap equation. We have shown, however, 
that this fact will not change much about our conclusions concerning the ghost propagator. 

Our major conclusion about the ghost propagator is obtained from the joint use of lattice 
data and of a result stemming from the ghost propagator Dyson-Schwinger equation: this 
result consists in a simple and cut-off independent formula for the ghost propagator depen- 
dence at small momentum. If we choose 1.5 GeV as the renormalisation scale, we get from 
lattice 

F (1.5 GeV) = Zg ~ 1.6 whence Fr (0, 1.5 GeV) ~ 2.2. (35) 

This has of course to be refined particularly regarding finite volume effects which should 
be considered with more care. The major point in this note, from the point of view of 
the renormalisability of the theory, is that lattice artefacts behave as powers of a (in this 
case 0(a^)). Of course at /3 = 5.7 lattice spacing is not yet small and this leads to a 
significant uncertainty which deserves further study. Any statement from lattice concerning 
bare quantities has to be taken with great caution since the very slow logarithmic dependence 
has chances to escape numerical observation. 
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